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Abstract. The effect of viscosity and of converging flows 
on the formation of blobs in the slow solar wind is analysed 
by means of resistive MHD simulations. The regions above 
coronal streamers where blobs are formed ( [Sheeley et aLj 
1997') are simulated using a model previously proposed by 
Einaudi et al.^^1999 ). The result of our investigation is two- 
fold. First, we demonstrate a new mechanism for enhanced 
momentum transfer between a forming blob and the fast so- 
lar wind surrounding it. The effect is caused by the longer 
range of the electric field caused by the tearing instability 
forming the blob. The electric field reaches into the fast so- 
lar wind and interacts with it, causing a viscous drag that is 
global in nature rather than local across fluid layers as it is 
the case in normal uncharged fluids (like water). Second, the 
presence of a magnetic cusp at the tip of a coronal helmet 
streamer causes a converging of the flows on the two sides of 
the streamer and a direct push of the forming island by the 
fast solar wind, resulting in a more efficient momentum ex- 
change. 



1 Introduction 

The understanding of slow solar wind genesis has progressed 
considerably in the last few years following the observational 
discovery of plasma inhomogeneities (called blobs) formed 



and expelled from regions above coronal streamers ( |Shee 
ley et al. 1997[ Wang et al. 1998| ). The observational ev 
idence has been provided by the Large Angle and Spectro 
metric Coronagraph (LASCO) instrument on the Solar and 
Heliospheric Observatory (SOHO). Approximately 4 blobs 
per day were observed in a relatively quiet period for the so- 
lar corona (year 1997, near solar minimum). After formation 
the blobs are carried away with the ambient slow solar wind. 
Generally, the plasma blobs move radially outward with their 
speed increasing from 0-250 km s~^ in the region 2-6 R© 
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(C2 field of view) to 200-450 km s"^ in the region 3.7-30 
(outer portion of the C3 field of view). 
Models of the formation of blobs have focused on mech- 
anisms involving reconnection between open field lines and 



closed field lines (Einaudi et al.| 1999; Einaudi et al. 2001 


Endeve et al.l ,2003 


Fisk and Schwadron, ,2001 


|van AalstI 


et al.| 119991 IWang et al.| |1998[ |Wu et al.| 2000 


). Through 



field line reconnection, plasma from the denser regions 
within the helmet streamer is liberated and becomes tied to 



open field lines and forms the slow solar wind ( Einaudi et al. 

[2Dog. 

The present work is particularly centred around the model 
proposed by [Einaudi et al.| ( p"999| ). In the region downstream 
of the cusp in the streamer belt, Einaudi et al. (1999) ap- 
proximate the configuration as a ID reversed field configu- 
ration with a ID wake velocity profile. Under these assump- 
tions, they can reproduce both the acceleration of the slow 
solar wind and the formation of blobs, formulating a com- 
plete picture of the genesis of the slow solar wind, in ac- 
cordance with the observations ( [Einaudi et al.[ |1999| ). The 



scenario has recently been further extended ( Rappazzo et al. 



2005| to include the so called melon seed effect due to the 



diamagnetic force caused by the overall magnetic field radial 



gradients in the spherical geometry of the Sun ( Schmidt and 



Cargill[[2000| ) and to include the effects of the cusp magnetic 



configurations atop the closed field lines in coronal stream- 



ers ( ^Lapenta and Knoll] |2005, ). 

The fundamental conclusion of the model by|Einaudi et al. 



(1999) and in successive refinements is that a blob is formed 
by the tearing instability and that by virtue of the interaction 
with the fast solar wind the blob is accelerated and ejected, 
carrying with it the plasma that forms the slow solar wind. 

The fundamental question posed by the present work is: 
"what is the basic physical process for the plasma accelera- 
tion"? We consider specifically two possibilities. First, that 
the momentum transfer causing the blob acceleration is due 
to viscous drag. Second, that the acceleration is due to the 
insurgence of the tearing instability and its non-linear ef- 
fects. The acute reader will not find the suspense spoiled 
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by learning that, as usual, Nature turns out not to be black 
or white. We discover that both mechanisms are present in 
a interplay mediated by the electric field. Furthermore, the 
presence of the converging flows at the cusp of the coronal 
streamer causes a direct momentum transfer by the action of 
the plasma flowing against the forming island and pushing it 
upward from the solar surface. 

The paper is organised as follows. Section 2 summarises 
the initial configurations considered distinguishing the ID 



2.1 ID force-free equilibrium 



initial configuration proposed by |Einaudi et al.| ( |1999| ) and 
its 2D extension to include the presence of a magnetic 
cusp (' Lapenta and Knoll 2005 1 ). Section 3 very briefly 
provides the methodology for our investigation: the re- 
sistive MHD model as implemented in the FLIP3D-MHD 
1991| ). Section 4 and 5 present the non- 



code (Brackbill 



linear evolution, respectively, of the case of the initial ID 
reversed field and of the 2D configuration with a magnetic 
cusp. Conclusions are reached in Section 6. 



2 Initial configuration 

The region where the blobs are formed and presumably re- 
connection happens has been localised thanks to the LASCO 
observations. The blobs are formed between 2-3 R© (the 
cusp region) and are accelerated afterwards, reaching termi- 
nal speed 250-450 km s~^ at 20-30 R©. The sonic point 
where the fast solar wind reaches sonic speeds is estimated 
at 5-6 R©, the Alfven point is further outwards. In the region 
of blob formation the fast wind is slightly subsonic and sub- 
alfvenic. In the present paper we consider the region where 
the plasma blobs and the slow solar wind originate. We con- 
sider the region near the cusp where we assume the plasma 
to be near sonic but subalfvenic. 

We consider two types of configurations previously pro- 
posed for the regions where blobs are observed to form. 
The first is ID and is a simple force-free extension of the 
Harris equilibrium ( [Einaudi et al.| |2001| ) 



The second, is 
a more complex 2D equilibrium representative of a Helmet 
streamer ( Lapenta and Knoll|,2005 ). In both cases, the cen- 
tral current sheet represents the area above the tip of a helmet 
streamer where the blobs have been observed to form in the 
LASCO images. To model the presence of the fast solar wind 
emanating from the coronal holes, the initial equilibrium in- 
cludes the presence of a flow along the field lines. The flow 
is chosen to be zero along the field lines directly above the 
streamer tip and to increase rapidly away from it, to reach an 
asymptotic value corresponding to the fast solar wind. The 
details of the initial configurations have been presented in 
previous publications and are summarised below for each of 
the two cases. 

We use Cartesian co-ordinates where y and z are parallel 
to the Sun's surface (with z along the meridians and y along 
the equator) and x is orthogonal to and pointed away from 
the Sun's surface. 



We use the ID model proposed by Einaudi et al. ( |2001 ). The 
magnetic field is force free (i.e. the pressure, temperature 
and density are uniform) and is given by 



Bqx = Bq tanh((z - zo)/L) 
Boy = Bosech{{z - zq)/L) 



(1) 



where L is the length scale of the variation of B. A plasma 
flow with a wake profile is chosen: 



uqj:{z) = uosech{{z - zq)/L); 



(2) 



The flow profile is chosen according to the model of [Ein- 



audi et alT] ( |2001| ) and represents the fast solar wind propa- 
gating along the open field lines, while a stationary plasma 
lays above the helmet streamer. As noted in Einaudi et al.J 
(2001) , the scale-length of variation of the flow shear might 
differ significantly from the length scale of variation of the 
magnetic field. However, for simplicity we assume them to 
be the same. In the present paper we choose the particular 
case uq/va = 2/3, where the Alfven speed va is defined 
with ^o- 

Periodic boundary conditions are applied in the x direc- 
tion and free- slip boundary conditions are applied in z. The 
boundaries in z allow no outflow of plasma and form ef- 
fectively a closed plasma channel. A small perturbation is 
added to the system to stimulate the growth of any instability 
present. The perturbation is added to the y component of the 
vector potential as: 

SAy = esm{kxx)sin{kzz) (3) 

with kx = 27r/Lx, k^ = tt/L^ and e = Bo/15. 

2.2 2D Helmet- streamer 

The second initial configuration is based on a recent work 



by Wiegelma nn et al. ( 1998) where a magnetic field config- 
uration with a current sheet typical of a helmet streamer is 
computed analytically. The equilibrium is rigourously exact 
only for infinitely stretched helmet streamers, but it can be 
used even for realistic aspect ratios as a good approximation. 

In the co-ordinate system used here, the current is initially 
directed along y and the initial magnetic field is in the (x, z) 
plane (see Fig. [T]), with x being the direction away from the 
Sun's surface. 

The initial equilibrium is 2D and is independent of the y 
co-ordinate. The initial configuration is characterised by the 
y component of the vector potential: 



Ay{x^ z) = — log < cosh 



Po 



{z - Zq) 



where c, zq and k are constants and Pq{x) = sie 
The initial magnetic configuration has a current: 



-log — 

c k 

(4) 

S3 



J 



cosh 



(5) 
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Fig. 1. Initial Helmet Streamer equilibrium. Note that the vertical 
and horizontal axis are not on the same scale, the vertical axis in 
compressed. 



and is held in pressure equilibrium by a pressure profile: 

Po 



P ■ 



cosh^ {^/fc{z - Zo)) 



Pb 



(6) 



The separatrices between open and close field lines are posi- 
tioned symmetrically at equal distances from the centre zq : 




arctanh 



Po 



Po 



(7) 



where = ke . In the simulations shown below we 
choose the parameters as ( [Wiegelmann et al.] |2QQ0| |1998[ 
Lapenta and Knoll| [2003] ) si .8, 52 4, 53 = .2, c = 15, 
As = .1073 and k = 1, zq = 1. Following [Wiegelmann 



et al. (2000), a small (20%) background plasma pressure pb 
is added to simplify the numerical treatment. The initial con- 
figuration is assumed to have uniform plasma density and the 
specific internal energy profile can be derived from eq.([6]): 



P 



P(7 - 1) 



Additionally, we add an initial wake velocity field: 



u{z) 



1 — sech 



\Z - Zo\ - Zs 

Ly 



(8) 



(9) 



where b is the unit vector in the direction of the magnetic 
field, where Ly is the scale-size of variation of the velocity 
field. The flow is present only outside of the two separatrices, 
where the field lines are open and is zero inside, where the 
field lines are closed. 

At the bottom boundary (x = 0), the field lines are tied and 
the flow and the plasma are continuously resupplied with the 
same value used for the initial condition to simulate the fast 
wind ejected from the coronal holes. At all other boundaries, 
open conditions are applied to allow the outflow of plasma. 



Unlike the previous case, no initial perturbation is applied 
since the error in the approximation of infinite stretching typ- 
ical of the derivations leading to eq. (|4]) is in itself a sufficient 
initial perturbation. 

In the present 2D equilibrium the equilibrium vary along 
z and the typical thickness of the layer downstream of the 
cusp, L = Lz /20. As in the ID case, the velocity is assumed 
to have the same length scale as the magnetic field: Ly = L. 



3 Simulation approach 

The simulations for the present work are conducted with the 
FLIP3D-MHD code ("BrackbiU' TMT) based on the viscous- 
resistive MHD model (Goossens , 2003^), comprising a mass 
continuity equation. 



dp 
'dt 



pV • u = 



Faraday's and Ampere's laws, 
B 



dt 



B 1 

= — • Vu - - [V X 7^ J] 

P P 



(10) 



(11) 



IIqJ = V X B 



Ohm's law. 



E = 7^ J - u X B 
a momentum equation. 



(12) 



P 



du. 

dt 



-V 



P- 



B^ 

2/io 



BB 

Po , 



-VApV-u- v-i^pn 

(13) 



and an energy equation. 



p- = -pV-u- 



- Ap(V • u)2 + iyp{U ■ n) + r]{J ■ J) (14) 



where p is the mass density, B is the magnetic field inten- 
sity, J is the current density, c is the speed of light, po is 
the magnetic permeability of vacuum, u is the fluid velocity, 
/ is the specific internal energy, and p is the fluid pressure. 
The symmetric rate-of- strain tensor, 11, is defined in the usual 



way (BraginskiH[1965[ ), 



n = -[Vu + Vu^] 



(15) 



The pressure is given by the equation of state, p = (7—1)/)/ 
with 7 = 5/3. The transport coefficients are the kinematic 
shear viscosity the kinematic bulk viscosity A, and the re- 
sistivity T]. The solenoidal condition on B, 



V-B = 



(16) 



is imposed throughout the evolution. 

The system under investigation is studied in a two- 
dimensional plane (x, z) using FLIP-MHD, a resistive MHD 



code described injBrackbill (1991 ). The system size in each 
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direction is labelled Lx, Lz and a grid made of 300 x 600 La- 
grangian markers is used in each direction arrayed initially 
in a 3 X 3 uniform formation in each cell of a 100 x 200 
grid. The units of the simulations are normalised to the mag- 
netic (and velocity ) scale length L and to the Alfven time 
ta = L/va- The Lundquist number ( [Goedbloed and Poedts 



2004 ) is computed a.s S = /hoVaL/t] and the Reynolds num- 
ber as: Ry = uopL/v (using the initial uniform density and 
the fast solar wind speed). We consider here for simplicity 
only shear viscosity (A = 0). 



4 Non-linear evolution of the ID initial equilibrium 

The initial state is characterised by both a flow shear and a 
magnetic shear. In principle, a number of instabilities could 
be present ( |Biskamp[ri997t|Bettarini et al.||2QQ6| ). However, 
for the choice of parameters considered only two types of 
instabilities are present. 

As shown in previously published results ( [Einaudi et al.[ 
|2QQ1[ [Lapenta and Knoll] |20Q5| ), the resistive tearing insta- 
bility is present and grows at a rate determined by the re- 
sistivity of the system. The resistive tearing mode leads to 
the formation of a magnetic island that mixes outer and inner 
field lines. In the process, the outer solar wind can transfer 
its momentum to the inner plasma and produce a net motion 
of the tearing island. This process is clearly reminiscent of 
the observed formation of blobs and their subsequent accel- 
eration away from the Sun. 

However, a second instability is present. Even in absence 
of resistive MHD modes, the viscosity is itself capable of al- 
tering the initial profile. Indeed, a viscous flow would tend to 
flatten any flow profile in time by virtue of direct viscous drag 
alone. The inner layer that is initially motionless would pro- 
gressively acquire speed ( jLandau and Lifshitz| |1959 ). In the 
corona, the viscosity is very low unless anomalous effects be- 
come important. However, in any realistically feasible sim- 
ulation, numerical viscosity has to be present to avoid the 
formation of small scale structures that cannot be resolved 
and would lead to numerical instabilities. 

One aim of the present work is to elucidate the relative role 
of this second process relative to the first. The question is: is 
the slow solar wind acceleration observed in the simulations 
due to the tearing instability or to the viscous drag? And how 
does it is scale with viscosity and resistivity, 

4. 1 Formation of the slow solar wind and of plasma blobs 

To exemplify the typical evolution observed in the simula- 
tions of the ID initial state, we report the results of a simula- 
tion for the particular case with Lundquist number = 10^ 
and viscous Reynolds number Ry = 10^. 

The initial phase of the evolution sees the formation of a 
magnetic island that perturbs the initial topology of the sys- 
tem and the initial flow profile. Figure [2] shows the velocity 
away from the Sun in false colour, superimposed on the mag- 
netic filed lines, at four times showing different phases of is- 



land formation. The island grows and moves away from its 
original formation site. Comparing the different times shown 
in Fig. |2] demonstrates that the island moves along with the 
flow. Note that in the simulations, the frame of reference is 
chosen stationary with the plasma flow in the flanks, so the 
motion of the island seen in the figure is with respect to the 
fast solar wind. 

The flow profile at different times is shown in Fig. [3] 
The momentum transfer between the flanks and the centre 
is made evident by the progressive flattening of the profile. 
The central plasma is progressively accelerated and the speed 
differential between the centre and the flanks diminishes with 
time. The temporal evolution of the speed of the central re- 
gion of the plasma is shown in Fig.|4]-a. Two cases are shown, 
besides the simulation considered above, another with higher 
viscosity (Ry = 10^) but equal resistivity is also reported. 
The process of acceleration progresses in two steps. After a 
progressive and continuous phase of acceleration of the ini- 
tially stagnant plasma, a sudden acceleration develops when 
the island becomes so big as to feel the effect of the boundary 
conditions. AS shown below, this last phase is only present in 
the low resistivity cases and is it is due to a transition from a 
X-point reconnection process to a Y-point reconnection pro- 
cess with a progressively more elongated leg'Biskamp|( |2Q00 



becoming progressively more turbulent and undergoing sec- 
ondary island formation. These features of the reconnection 
mechanism are well known and will not be addressed here 
again, the interested reader is referred to the excellent text- 
book by [Biskampl ( [2000 ) 



The process of island formation is further shown in Fig.|4]- 
b where for the same two simulations, the reconnected flux 
is shown. The onset of the faster reconnection phase also is 
similarly due to the boundary conditions chocking the flow 
which turns inward toward the reconnection site at the x- 
point formed by the island driving faster reconnection. Again 
this process is an artefact of the boundary conditions. 

Figure |4] shows a comparison of two cases with different 
viscosity and demonstrates that indeed viscosity has a strong 
impact on both island formation and slow solar wind accel- 
eration. In the next paragraph this effect is further analysed. 

4.2 Role of Viscosity and Resistivity 

The behaviour displayed above is in agreement with previ- 
ously published results ( [Einaudi et al.| [200 1[ [Lapenta and 
Knoll 2005 ) and it is reproduced primarily to provide the 
reader with a summary of the typical evolution. Below we 
investigate the main point of the present effort, namely what 
are the causes of the observed behaviour displayed in Fig. [3] 
Both the tearing instability and viscous stresses can be the 
cause: to determine which dominates, we change indepen- 
dently viscosity and resistivity. The tearing instability is pro- 
portional to resistivity, viscous drag is proportional to viscos- 
ity. We need to determine whether the momentum transfer is 
proportional to viscosity or resistivity. 

Figure [5] shows the summary of the different runs con- 
ducted varying resistivity and viscosity. Each panel corre- 



G. Lapenta and A.L. Restante: Blob formation in the solar wind 



5 




sponds to a different resistivity and reports the time evo- 
lution of the maximum difference in the averaged profile 
< Ux{z) >x, corresponding to the differential between the 
flank velocity (fast solar wind) and the central plasma (slow 
solar wind). In each panel several runs are displayed, corre- 
sponding to different viscosities in each run. 

Two trends are clear. First, for some runs, the same pro- 
gressive acceleration seen in Fig. |4] above is followed by a 
disruptive phase when the boundary conditions affirm their 
presence. In some cases at low Ry, the system never reaches 
this state before the completion of the run. Second, observing 
the progressive acceleration phase, one unmistakable conclu- 
sion emerges: the speed of acceleration (i.e. the slope of the 
curve in the figures) is hardly affected by the resistivity but 
scales monotonically and markedly with the viscosity. The 
conclusion regarding the comparatively lower sensitivity of 
the acceleration mechanism to resistivity was already dis- 
cussed in a previous work ( Lapenta and Knoll) |2QQ5| ). Here 
the focus is on the much more apparent sensitivity to vis- 
cosity. As viscosity is increased by two orders of magnitude 



from a Reynolds number of Ry = 10"^ to Ry = 10^, the 
slope increases markedly: for example, the speed at time of 
t/TA = 200 for the runs with S = 10^ drops by 10% in the 
Ry = 10^ case and 55% in the Ry = 10^ case. Clearly the 
mechanism for slow wind acceleration appears to be more 
linked to the process of viscous drag than to the processes 
related to resistivity and island formation. 

The either/or test just conducted seemingly appears to con- 
clude that the main mechanism in action in the simulations 
is just simple viscous drag. However, a closer scrutiny re- 
veals that this conclusion based on the overall evolution of 
the systems still misses one important contribution that the 
island formation provides. The next subsection reveals this 
more subtle effect. 

4.3 Role of the Electric Field 

A key aspect of the evolution of the system considered above 
rests in the role of the electric field. The initial equilibrium 
is already characterised by an electrostatic field needed to 
produce the velocity shear imposed initially. The subsequent 
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A) 




B) 



Re=10^ 



50 100 150 200 250 300 350 400 450 



Fig. 3. Evolution of the average velocity profile < Vx{z) >x (nor- 
malised to the asymptotic speed uq) for two runs starting from the 
initial force free reversed field equilibrium with the same resistivity 
(S = 10^) but different viscosity (panel A with Rv = 10^, panel B 
= 10^). 

evolution of the tearing mode leads to the growth also of an 
electromagnetic field. 

In the present MHD treatment, the plasma flow is primar- 
ily due to the drift motion caused by the presence of electric 
fields: 

ExB 

VExB = (17) 

Once the electric field is written in terms of the vector and 
scalar potential, the drift clearly shows the presence of a dual 
nature: 

1 dA 

VEXB = -^(^^ ''^^^''^^ ^^^^ 

The flow can be of two natures: electrostatic and electromag- 
netic. The initial configuration has only 8iVx{z) component 
of the flow that is caused by a component of the electric 
field: therefore the electric field has zero curl and it is purely 
electrostatic. The tearing mode is a mode primarily involv- 
ing the evolution of the out of plane component of the vector 



Fig. 4. Evolution of the peak of the average velocity (normalized 
to the asymptotic speed m) max2(< Ux(z) >x) (A) and of the 
reconnected flux (B) for two runs starting from the initial force free 
reversed field equilibrium with the same resistivity (S = 10^) but 
different viscosity (panel A with Ry = 10^, panel B Rv = 10^). 

potential A and it is inductive in nature since it produces an 
electromagnetic field thanks to the topological variations of 
B. However, in presence of out of plane components of the 
field and of flow shears, there is also an electrostatic potential 
associate with it ( [Daughton and Karimabadi[ |2QQ5 ). 

The viscous effect alone, instead, does not directly affect 
the magnetic field and causes electrostatic fields. If we con- 
sider the evolution of the initial state under the action of vis- 
cosity alone, the full set of equations reduces to just: 

dux 1 d'^Ux 

— - = (19) 

dt Ry dz 

governing the only remaining variable Ux{z). Using Ohm's 
law: 

V X £; = -VV • B + BV • V - B • Vv + V • VB (20) 

where all terms on the right-hand side are zero: the first be- 
cause of the absence of magnetic mono-poles, the second be- 
cause the only non-zero component is Ux which is only a 



G. Lapenta and A.L. Restante: Blob formation in the solar wind 

A) 5' = 10^ B)S = 2- 10^ 






Fig. 5. Evolution of the difference AVx between the maximum 
(fast solar wind) and minimum (slow solar wind) in the profile of 
the average vertical velocity (normalized to the asymptotic speed 
uo) ma^z{< Ux(z) >x)- Several runs are shown starting from 
the same initial force free reversed field equilibrium, each panel 
corresponds to a different resistivity: A) S ^ 10^, B) 5* = 2 • 10^, 
C) S = 10^ and D) 5* = 2 • 10^. In each panel 5 viscosities are 
shown, corresponding in order from the top curve to the bottom 
curve to Ry 10^, Ry = 2 - 10^ Ry 10^ Ry = 2 ■ 10^ Ry 
10^. Note that in the legends, the viscosity is reported directly rather 
than the Reynolds number: in our normalized units they are simply 
one the reciprocal of the other. 



function of z (and therefore the divergence of u is zero), the 
third because is constant along field lines, and the last 
because, similarly, B does not change along the flow lines. 
The variation of the velocity field generated by viscous drag 
causes an electric field with zero curl, i.e. it generates an 
electrostatic field. 

The effect of viscosity alone can be computed easily by 
solving directly eq. (19). As it is well known ( [Barenblatt| 
,1996^ ), a linear diffusion equation like eq. (19) admits a self- 
similar analytical solution of the form: 



1 



(21) 



that provides a simple scaling law for the time-variation of 
the differential in the speed of the central initially stagnant 
plasma and the outer fast solar wind: AVx ^ 1/ V^- 

To illustrate the effect of viscosity alone, we have solved 



the equation for viscous drag (i.e. eq. ( 19 )) numerically start- 
ing with the same initial profile considered in the full MHD 
simulations and for the same ranges of viscosity. Figure |6] 
shows the evolution in time of the jump in speed between the 
central plasma and the outer fast solar wind. The results can 
be compared directly with Fig. |5] for the full MHD simula- 
tions. As can be observed the additional effect of the pres- 
ence of the magnetic field evolution is very significant. The 
general shape of the time evolution is the same in the purely 
viscous case as in the full MHD case, and in both cases re- 




Fig. 6. Evolution of the difference in velocity AVx between the 
central plasma (initially at rest) and the outer fast solar wind. The 
results are from a simple purely viscous simulation that neglects any 
magnetic field and resistive effect. Five viscosities are shown, from 
top to bottom: Ry = 10^, Ry = 5 ■ 10^ Ry = 10^ Ry = 5 ■ 10^ 
and Ry = 10^ 



sembles the inverse-square-root power law predicted by the 
self-similar solution. A quantitative comparison reveals that 
the evolution with full MHD has a very much increased rate 
of decay. For example, considering the case Ry = 10^, at 
time t/rA = 200, the full MHD simulation shows a differ- 
ence in speed of 0.65, while in the purely viscous case this 
difference is still 0.87. Clearly, the full MHD evolution pro- 
vides additional drag forces that transcend the direct viscous 
drag. 

The additional force is coming from the electric field. 

Figure [7] shows the vector and scalar potential at time 
I/ta = 76, for the same simulation shown in Figs. |2] above. 
A significant electrostatic field is present form the beginning 
and it is caused, as noted above, by the need to support the 
initial sheared velocity field. The history of the electric field 
energy is shown in Fig. [8] The partial contribution of the 
electrostatic field is shown as a dashed line. Except for the z 
component that includes the initial strong electrostatic field 
supporting the sheared flow, the electromagnetic component 
proper of the tearing mode dominates all components. The 
electric field perturbation is therefore primarily caused by the 
tearing growth and not by viscous drag. 

The tearing instability developing in the centre of the sys- 
tem moves with respect to the plasma in the flanks as shown 
in Fig. [2] The island remain localised in the centre of the 
plasma and it comoves with the local plasma speed. The di- 
rect viscous drag is small. However the island's electric field 
structure extends far beyond it and well into the flanks, as 
shown in Fig. [7] These fields move with the island speed and 
have a relative speed with respect to the fast solar wind in the 
flanks. The expanded range of influence of the island deter- 
mined by the electric fields allows for an enhanced drag that 
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Fig. 7. Vector and scalar potentials at time t = 76 for a run starting 
from the initial force free reversed field equilibrium with 5* = 10^ 
and Rv — 10^ (same run considered in Fig. [5]). 




Fig. 8. History of the electric field energy for a run with the ini- 
tial force free reversed field equilibrium and with 5* = 10^ and 
Rv — 10^. The electrostatic component is shown as a dashed line, 
the total as a solid line. Note that in a 2D simulation in the (x, z) 
plane only the two in-plane components can be present for the elec- 
trostatic field. 



no longer happens only across layers touching each other but 
extends globally as the electric field of the island feels the 
drag force over the whole domain. 



5 Non-linear evolution of the initial 2D helmet streamer 
equilibrium 

The 2D helmet streamer configuration considered is charac- 
terised by the presence of a cusp region above the helmet 
streamer where the field lines from converging towards the 
cusp at lower altitudes become essentially parallel. In this 
region the field-aligned speed also changes direction from 



a converging motion at lower altitudes to a parallel motion 
above. In the ID equilibrium considered above, the focus 
was entirely on the region above the cusp where the flow and 
filed lines are initially straight and parallel. The 2D equilib- 
rium with cusp, instead, allows us to consider also the motion 
in the lower part where the flow and the field lines are con- 
verging. 

Converging motions are particularly relevant to the forma- 



tion of blobs observed in the LASCO images ( Wang et al. 



1998| : converging motions cause magnetic fluxes to pile 



up towards the converging point and can drive magnetic re- 
connection. Magnetic reconnection driven by macroscopic 
flows becomes enslaved to the flows and its rate no longer 
depends on the details of the microscopic physics allowing 
reconnection (Biskampl[200Q||Dorelli and Bim| [20031 [KnoIT 



and Chacon , 2002 k A previous study showed that in the 2D 



configuration under consideration the flow towards the cusp 
produces the conditions for a reconnection process that pro- 
gresses at the speed dictated by the incoming flux, insensitive 
to the resistive processes allowing the break up of the frozen- 
in condition ([Lapenta and KnolI| [2005|). As illustrated by 



[Lapenta and 



Cnoll ( 2005] ), the mechanism allowing the de- 



coupling of resistivity and reconnection rate is the flux pile 
up: flux piles up in front of the diffusion region to bring the 
speed of reconnection up to the needed rate ( Dorelli and Birn, 
2003] ). Indeed, increasing the magnetic field strength at the 



inflow side of the diffusion region increases the Alfven speed 
and speeds up the whole process of reconnection. 

In the present study, we consider, instead, the role of vis- 
cosity. As the brief summary of the previous understand- 
ing implies, the expected outcome is that in driven reconnec- 
tion viscosity should be a diminishing factor: as viscosity 
increases, the flows that drive reconnection become impeded 
by the viscous drag and see their ability to drive reconnection 
diminished. Therefore, the expectation is that viscosity is not 
a needed ingredient in the physics leading to blob formation 
in the region around the cusp. 

To investigate this point. Fig. [9] shows the magnetic sur- 
faces after the formation of a magnetic island above the cusp 
of the helmet streamer. As can be observed, indeed the pro- 
cess happens regardless of viscosity, and, as expected, the 
more viscous run shows a diminished speed of formation of 
the island. At a given same time in the two runs, the island is 
less developed in the higher viscosity case. 



This conclusion is further confirmed by Fig. 10 that shows 
the evolution of the reconnected flux for the two runs with 
different viscosity (but equal resistivity) The reconnection 
process is similar but slightly faster for the lower viscosity 
run. 

However, the fundamental difference of the 2D equilib- 
rium when compared with the ID equilibrium of the previous 
section is that also the acceleration of the blob is independent 
of viscosity. The mechanism allowing the acceleration is no 
longer the viscous drag exerted by the outer plasma on the 
forming island (and mediated by the electric field) but rather 
is directly the converging flow at the cusp. 
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Fig. 9. Non-linear evolution of the initial 2D helmet streamer equi- 
librium. Contours of the vector potential Ay, representing magnetic 
surfaces, at time IIta — 200 



Figure 1 1 shows the vertical speed in the same two runs 
with different viscosity considered above. The flow velocity 
originating from the open field lines converges and focuses 
at the cusp in the region right below the forming island. The 
island finds itself in the same situation of a ball sitting atop 
a jet in a fountain: it is pushed by the jetting plasma. The 
converging flow focuses at the cusp and pushes the island 
upward by direct momentum transfer from the plasma to the 
island field lines. Viscous drag is no longer needed. 



This last conclusion is further confirmed by Fig. 12 that 
shows the acceleration of the initially stagnant plasma above 
the cusp as a function of time. Clearly, viscosity is not the 
cause of the acceleration but rather just as it slows down by 
a small but perceptible amount the process of reconnection, 
it also slows down the ability of the converging flow to focus 
at the cusp and push the forming island upward. 




Fig. 10. Non-linear evolution of the initial 2D helmet streamer equi- 
librium. Comparison of the time evolution for the reconnected flux 
in two runs with different viscosities {R^ — 10"^ and R^ — 10^) 
but equal resistivity S — 10^. 



6 Conclusions 

We have investigated the role of viscosity in the physical 
mechanisms responsible for the formation of blobs in the 
slow solar wind. We started from two types of models. The 



first is a initial ID equilibrium proposed by |Einaudi et al. 
( 2QQ1| ) and designed to capture the field reversal and the flow 
shear present in the plasma above the helmet streamers in 
the regions where blobs have been observed to form in the 
LASCO images ( [Wang etaL 1998 ). The second model is 
based on a 2D equilibrium that extends the previous case by 
including the effect of converging flows at the cusp imme- 
diately above a helmet streamer ( Lapenta and Knoll 2005 1 ). 
The role of viscosity is different in the two cases. 

In the ID initial equilibrium, viscosity is shown to provide 
the basic mechanism for momentum transfer from the slow 
solar wind to the forming blob. The viscous drag produces a 
progressive acceleration of the plasma blob. The rate of in- 
crease of the slow solar wind speed (acceleration) is propor- 
tional to the viscosity in the simulation. However, an unex- 
pected finding is that viscosity is not acting alone in its usual 
local drag action across layers of flowing plasma. Rather, the 
coupling between the accelerating island and the fast solar 
wind in the flanks is mediated by the global electric fields 
caused by the tearing instability forming the island. A direct 
numerical comparison of the evolution under pure viscous 
action and the full MHD evolution including the generation 
of electromagnetic fields due to the tearing instability proves 
that the drag action is much increased from simple viscous 
drag. 

In the 2D initial configuration, the momentum transfer ac- 
tion is completely different. The viscous drag is not a de- 
termining factor in that it does not provide the mechanism 
for momentum transfer from the fast wind to the forming is- 
land. In the 2D case, the flow merges at the cusp and simply 
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Fig. 11. Non-linear evolution of the initial 2D helmet streamer equi- 
librium. False colour representation of the vertical speed away from 
the Sun Ux, at time IIta — 200. Two runs are shown with different 



viscosities (panel K\Rv 
resistivity S — 10^. 



lO"" and panel B: Ry = 10^) but equal 



pushes the island by direct momentum transfer acting in the 
same way a ball is pushed up by jetting water in a fountain. 
The momentum transfer in this case is coming directly by 
the action of plasma particles directed toward the magnetic 
filed lines of the island and being deflected by them. Their 
change in momentum of deflected particles provides a direct 
accelerating force. Viscosity has no role to play and its effect 
is just to diminish the flows and to reduce the rate of island 
formation and its subsequent acceleration. 

The results presented are relevant to the Sun where viscos- 
ity and resistivity due to direct collisional transport are ex- 
tremely low. The ability to explain both in the ID and in the 
2D case how momentum transfer is provided by extra effects 
not directly due to viscosity is crucial. The results obtained 
here, as all other results obtainable within the range of cur- 
rently achievable simulation cannot properly use the viscos- 
ity and resistivity of the real corona and cannot account self- 
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Fig. 12. Comparison of the time evolution of the vertical velocity 
in the central axis in two runs starting from the initial 2D helmet 
streamer equilibrium for the maximum with different viscosities 
{R^ = 10^ and R^ = 10^) but equal resistivity S = 10^. 



consistently for anomalous processes. It is therefore crucial 
to show that the accelerating mechanisms presented here are 
indeed stronger than purely viscous processes would produce 
and can sustain the scaling to realistic coronal parameters. 

We remind the reader that we conducted the research here 
on a Cartesian geometry where the melon seed effect caused 
by the expanding geometry around the Sun is not consid- 
ered (Rappazzo et al.| 2005| ): the results of the present work 
need to be integrated by considering that additional effect. 
Future work will consider how the conclusions reached here 
are modified by additional geometric effects. 
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